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Impact problem for the quasi-linear viscoelastic 
standard solid model 
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Abstract: The one-dimensional impact problem in the case of Fung’s quasi-linear 
viscoelastic model is studied for the relaxation function of the standard solid model (or 
Zener model). At that, quasi-linear viscoelastic Maxwell and Kelvin-Voigt models are 
recovered as limit cases. The results of numerical simulations for some illustrative values 
of the dimensionless problem parameters are presented. 
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1 Introduction 


Impact testing of soft biomedical materials and biological tissues represents a considerable 
practical interest ( |Kalcioglu et ah , 2011). In recent years there has also been a growing 
interest in impact testing at nanoscale (Constantinides et ah, 2008). The mechanical 


behavior of articular cartilage as a biological tissue demonstrates several complex fea¬ 


tures including viscoelasticity and nonlinearity (Hayes et ah, 1971), and there have been 


developed a number of mathematical models for soft biological tissues and biomaterials. 


including poroelasticity (Mow et ah, 1980), poroviscoelasticity (Setton et ah, 1993), and 


rheological network modelling ( Bischoff et al.| 2004). Due to the short time span involved 
in the impact problem, disallowing macroscopic motion of the interstitial fluid relative 
to the solid matrix of the tissue, we ignore its multiphasic nature. However, it should 
be noted that the validity of the assumption of viscoelasticity may be jeopardized if the 
indenter has a submicron diameter, because the macroscopic size of interest has been 
reduced to a very short distance. 

Needless to say that the interpretation of experimental data on material characteri¬ 
zation by high-rate impact tests crucially depends on the mathematical material model 
employed for the analysis. Because biological materials exhibit time-dependent response 
to external mechanical stimuli, viscoelastic models seem to be appropriate candidates. 
The corresponding impact problem in spite of being simply formulated imposes consid¬ 
erable mathematical difficulties for analytical solution. The analytical solution of one¬ 


dimensional impact problem is known for the Kelvin-Voigt (Wineman and Rajagopal 
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2000 Popov, [2010 ) and Maxwell (Stronge, 2000) models, and in the case of viscoelas¬ 


tic standard solid model it was studied by means of asymptotic methods (Butcher and 


Segalman, 2000; Argatov, 2013). Furthermore, some qualitative properties of the one¬ 


dimensional impact were established (Argat^ 2013) in the framework of general linear 
viscoelastic model represented by the Boltzmann’s hereditary integral. 

It is clear that linear impact models are not capable of describing sensitive features of 
impact phenomena for soft tissues (Edelsten et ah, 2010 Kalcioglu et ah, 2011). One of 


the most striking examples is the behavior of the coefficient of restitution as a function 


of the velocity of impactor observed in recent experiments on articular cartilage (Varga 


et al., 2007 Edelsten et ah, 2010), whereas the linear viscoelastic impact models asserts 


that it should be independent of the impact velocity. That is why, it was hypothesized 
(Edelsten et al., 2010 Argatov, 2013| ) that utilizing non-linear viscoelastic models in the 
impact problem may shed light on the behavior of the output impact variables (e.g., 
coefficient of restitution) with variation of the input impact parameters (impactor mass 
and velocity). 

In the present paper, we employ the quasi-linear viscoelastic (QLV) model, which was 
proposed by Fung (1981) and is based on the hereditary integral with a certain relaxation 
function. In particular, we consider the following relaxation functions: 1) Standard solid 
model (Zener model), 2) Maxwell model, 3) Kelvin-Voigt model. At that, the QLV stan¬ 
dard solid model is characterized by the following four parameters: Instantaneous elastic 
modulus, Eq, characteristic relaxation time, r^, equilibrium-to-instantaneous modulus ra¬ 
tio, p, and (dimensionless) nonlinearity parameter, B. Correspondingly, the QLV Maxwell 
and Kelvin-Voigt models are obtained from the QLV Zener model as limit cases as p —)■ 0 
(after appropriate renormalization). On the other hand, the QLV Fung model is an ex¬ 
tension of the Boltzmann integral to account for elastic non-linearity via the parameter 
B and therefore by dehnition the QLV Zener model reduces to a standard solid model as 
B ^0. 

The paper is organized as follows. In Sectioi][^ we give the general background for the 
QLV standard solid model based one-dimensional impact problem. To facilitate numerical 
implementation, the impact problem is reduced a system of the hrst-order differential 
equations. The two limit cases for the QLV Zener model are considered in Section]^ The 
results of numerical simulations for some illustrative values of the dimensionless problem 
parameters are presented and discussed in Sections and respectively. 
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2 Impact problem formulation 


2.1 Equation of motion of the rigid impactor 


We consider the schematic representation for impact loading of a biological tissne specimen 
shown in Fig. 2d^ On the basis of Newton’s second law, the impactor motion is determined 
by the differential equation 

mx =—F, (1) 

with the initial conditions 

a;(0) = 0, i;(0) = Uq. (2) 


Here, m is the mass of the impactor, F is the reaction force of the tissue specimen 
(which is assumed to be positive), Vq is the impactor speed at incidence, tc is the contact 
duration. In other words, tc denotes the instant when the specimen reaction force changes 
its sign and the impactor acceleration, x, vanishes, so that 


F 


t = tc 


= 0 . 


( 3 ) 
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Figure 1: Scheme of impact loading of the tissue specimen: (a) Initial conhguration before 
the contact; (b) Current conhguration. 


An important characteristics of the impact problem ([^, ([^ is the coefficient of resti¬ 
tution, e*, which is identihed as the ratio of the absolute value of the impactor velocity 
at separation, |h(tc)|, to the impactor speed at the incidence, i:(0), that is 


i{tc)\ 

Vo 


( 4 ) 


Moreover, some other impact characteristics also represent a considerable practical 
interest. Namely, the peak value, Xm, of the impactor displacement, which occurs at the 
time moment t = tm when x{tm) = 0, and the peak value, of the contact force, which 
is achieved at the instant t = Im, when F{tM) = 0. 
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2.2 Quasi-linear viscoelastic standard solid model 


We represent the specimen reaction force based on the qnasi-linear viscoelastic (QLV) 
model introdnced by Fnng (1981). The QLV model was proposed for modeling biological 
tissnes in one-dimensional stress state, which seems to be appropriate for blnnt impact 


testing (Edelsten et ah, 2010 Varga et al., 2007). Let ns recall that the one-dimensional 


stress, in the impacted specimen is determined as a ratio between the applied force 
— F{t) and the nndeformed cross-sectional area A. The corresponding extensional strain, 
e(t), is dehned as 

£(() = ^ - 1. (5) 


where h{t) is the thickness of the tissue specimen at the current time, and ho = h(0). 

According to the Fung QLV model, the dependence of the stress a{t) at time t on the 
strain history is given by 


t 

/X f X de , . , 

a(t)^ j K(t-s)—-(s)ds, 

0 - 


( 6 ) 


where it is assumed that e{t) = 0, t G (—cxd, 0). The lower integration limit 0“ indicates 
that the integration in (|^ starts at inhnitesimally negative time so that to include the 
strain discontinuity at time zero. 


Based on Fung’s formula (Fungi 1981 Rajagopal et ah, 2007 Muliana and Rajagopal 


2012) for a nonlinear elastic stress function a^[e{t)), the derivative da^/de, appearing in 


Eq. (^ is evaluated as 


^ = Eoexp(R£), 
de 


(7) 


where Eq and B are material parameters. Note that the elastic modulus Eq is related to 
the instantaneous elastic response of the tissue material. 

Correspondingly, the stress-relaxation function, K(t), which enters Eq. ([^ will be 
assumed to be normalized as 

77 ( 0 ) = 1 . ( 8 ) 


Due to geometric compatibility, the current thickness of the specimen, h{t), and the 
initial thickness of the specimen, ho, both entering Eq. ([^, are related to the impactor 
displacement, x{t), through 

Hq = h{t) + x{t). (9) 
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Therefore, in view of ([^, 0, and ([^, Eq. 0 yields 




Bx{s) \ dx ^ 
0 


( 10 ) 


In the present paper, the normalized stress-relaxation function is taken according to 


the standard solid model, also known as the Zener model, as follows (Wineman and 


Rajagopal 


Kit) = p + (I- p) exp ( - —) 

V TrJ 


( 11 ) 


Here, p represents the ratio between long-term and instantaneous responses, tr is a char¬ 
acteristic relaxation time. 

2.3 Reduction of the impact problem to a system of the first- 
order differential equations 

In view of (10), Eq. ([^ can be recast in the form 


mxit) = -^ I K{t-s)exp(^^)x{s)ds. 
ho J \ ho y 


Bxis) 


( 12 ) 


Now, by introducing non-dimensional variables 

t ^ X 

T = —, f =-. 

Tr VqTr' 


(13) 


the equation of motion (12) and the initial conditions (|^ can be rewritten as follows: 

+ aB = 0, (14) 

e(0) = 0, r(0) = l. (15) 

Here, the differentiation with respect to the dimensionless time variable r is denoted by 
stroke, B is the nondimensionalized reaction force given by 


f l^i'r - s)exp{(3^{s))^{s)ds, 


(16) 
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( 17 ) 


C(t) = p + (1 - p)exp(-r), 
and we introduced the notation 

AEqTr^ BvqTr 

a= - , —, 13 = —. —. 

TTIiIq IIq 

Furthermore, let us introduce an auxiliary notation 

T 

/(?.’■)= f exp(P({s))'^(s) ds, 


(18) 


0 - 


e r = 


dr 


(19) 


( 20 ) 


To proceed, hrst, we rewrite Eq. (16) with (17), (19), and (20) taken into account as 
follows: 


■7^(?.t) = p/(?,r) + (1 - rtexp(-r) / exp(s)exp(^5(s))e(s)<is. 


( 21 ) 


Now, differentiating Eqs. (19) and (21) with respect to r, we get 

= e(^)exp(/?^(r)), 

= e(r) exp(/?^(r)) - r) + p/(^, r). 

So, taking into account the above relations, we reduce the second order differential 
equation (14) to the following system of hrst-order differential equations: 

= 6, 

e' = —aE, 


< 


7' = eexp(/3^) - 7 + pf, 

J' = eexp(/?0- 


( 22 ) 


The initial conditions (15) are supplemented with two new zero initial conditions for 


the new variables 7 and /, that is 

«0) = 0, e(0) = l, ^(0) = 0, /(0)=0. 


(23) 


At that, in the non-dimensional variables (13), Eqs. ([^ and (|^ for the contact dura¬ 
tion, tc = trTc, and the coefficient restitution, respectively, take the form 


3F\ = 0 , 

'r=T-, ' 


( 24 ) 
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e. = -e(Tc)- 


(25) 


Note that the other impact characteristics introduced in Section 2.1 can be evaluated 


as 


Cm — tr'Tm, -Cm — ^I 


ha 


where Tm and tm are roots of the equations 


e|r=r„ = 0, J^'\ = 0. 

' ^ ’ I r=TM 


(26) 

(27) 

(28) 


Remark 1. The obtained system (22) contains two dimensionless parameters a and 
(3. In particular, the value of the parameter [3 determines the extent of nonlinearity. So, 

(29) 


if /3 = 0, then the Cauchy problem (22), (23) reduces to the system 

= e, e' = -a J", = e - T + p^, 

e(0) = 0, 6(0) = 1, .F(0)=0, 

which, in turn, results in the following differential equation of the third order: 

f!" + + ct^' + expf = 0, (30) 

subjected to the initial conditions 

e(0) = 0, e'(0) = l, e"(0) = 0. (31) 

The obtained Cauchy problem describes the linear impact problem for the Zener model 
studied in detail by Butcher and Segalman (2000); Argatov (2013). 


3 Limit cases for the QLV standard solid model 

3.1 QLV Maxwell model 

By passing to the limit as p —0 in formula ( pN| for the stress-relaxation function of the 
Zener model, we obtain the relaxation function for viscoelastic Maxwell model, Kuif), 
which in the non-dimensional variables ([T^ takes the form 

/Cm(^) = exp(-r). (32) 
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Correspondingly, formula (10) yields the following expression for the specimen reaction 
force: 

where 

ko = 


AEo 


h(] 


dE 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


In view of ([^ and^^, the non-dimensional contact duration, Tc, and the coefficient 
restitution, e*, are determined by the equations 

-T'mL =0, e* = -e(rc). (39) 


^m(C,^) = exp(-r) j exp(s)exp(/?^(s)) —(s)cis. 

0 - 

Similarly to Section |2.3| , it can be shown that the impact problem 

C + aEu = 0, m=E r(0) = l 

is equivalent to the system 

'e' = C 

= eexp(/3^) - Em 

with the initial conditions 

e(0) = 0, e(0) = l, -Fm(0) = 0. 

Note that the system (37) can be obtained from (|2^ by passing to the limit p = 0. 


Remark 2. In the limit case /3 = 0, the impact problem (31), (38) reduces to the 
Cauchy problem 

= e, e' = -aEu, E'm = e - Em, (40) 

e(0)=0, r;(0) = l, -Fm(0) = 0. (41) 


In this case, the coefficient restitution is given by the following formula (Butcher and 


Segalman, 2000; Strange. 2000): 


e* = exp 


-TT 


,_,. (42) 

— id 

Note also that the guantity l/{2ffa) has a physical meaning of the loss factor in the 
Maxwell model, which is schematically modeled as the combination of a spring (with the 
stiffness ko given by (34)) and a dashpot (with the damping constant b = r^ko) in serial 
connection. 



















3.2 QLV Kelvin—Voigt model 

First of all, let us rewrite formula ( |IT| in the form 

EoK{t) = +^—^exp (^- (43) 


where the elastic modulus E^o is related to the equilibrium elastic response of the tissue 
material. 

In the framework of the standard solid model, the following relation holds between 
Eoc and Eq: 

^oo P^O- 

Let us also introduce the so-called retardation time by 


P ' 


= 


Therefore, that Eqs. (10) and (0 can be recast as follows: 






X'(t) = l + i—^exp(-^). 

P \ PV 

Now, passing to the limit as p —?■ 0, the above formulas lead to the quasi-linear Kelvin- 
Voigt model: 

EKv{x,t) = ^^^ j iFKv(t- s)exp (44) 

0 - 


EKv{t) = 1 + 

Here, 6{x) is the Dirac delta function. 


(45) 


Thus, by substituting (45) into Eq. (|44|), we arrive at the relation 

ho 


EKY{x,t) = 


AE, 

ho 


(Exit) 

B Dp l-feT 




(46) 
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Further, introducing the non-dimensional variables 


z = exp 


I hn 


t 

1, '?■ = —> 
'^R 


(47) 


we 


reduce the impact problem Q, ([^, (46) to the following Cauchy problem: 

~/2 


Z"- 


= -a\z + z'){z + 1), 


^ + 1 

^(0) = 0, z'{0) = P'. 


Here we introduced the notation 


a = 


mhn 


/9' = 


BvqTr 

hn 


(48) 

(49) 

(50) 


According to (|^ and(|^, the contact contact duration, tc = t^jTc, and the coefficient 
restitution are identihed as 

^ + ^'L= =0, (51) 


e* — 


(3'z + l 


(52) 


Finally, the peak value of the impactor displacement, Xm, which occurs at the time 
moment tm = and the peak value of the contact force, Fm, which occurs at the 

time moment tu = Tr^m, can be evaluated as follows: 


h 

^'\r=r^ = 0 ’ Xm = + 1 ), 


z' + — 0, Fm — 

I t=tm 


AE, 

B 


(z(tm) + Z{tm)). 


(53) 


(54) 


Remark 3. Note that the change of the variable (47)i does not allow us to pass to 
the limit fd' = t) in the impact problem (51), (52). However, in the limit case (3' = D 
(that is when B = D), the impact problem under consideration reduces to the following 
Cauchy problem (see formula (46)), which corresponds to the linear viscoelastic Kelvin- 
Voigt model: 

A.E 

mx = - -^{x{t) + T'j^x{t)], te{0,tc), 

IIq 

a;(0) = 0, x(0) = uq- 
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In this case, the coefficient restitution is given by the following formula (Butcher and 


Segalman, 2000; Popov, 2010; Wineman and Rajagopal, 2000): 


e* = exp 




atari 


\/4 — 


a' 


a' 


(65) 


where \/~a'/2 is the loss factor in the Kelvin-Voigt model, which is schematically modeled 
as a parallel combination of a linear spring with the stiffness k^o = AE^o/ho and a dashpot 
with the damping coefficient b = ko 


''oqTr- 


4 Behavior of the main impact parameters 

First of all, to illustrate numerical solutions of the impact problems for the QLV Maxwell 
and Kelvin-Voigt models, we consider the following ranges for the model dimensionless 
parameters: f and f from 0 to 15, a from 0.25 to 100, and a' from 0 to 4. Note that these 


parameters, which are introduced in (18) and (50), can be viewed as non-dimensionalized 


reciprocal mass of the impactor and its initial velocity. 

In all hgures below, the following convention is used: (0) dot-dash black line represents 
the limiting case for the corresponding second parameter, while the solid lines (1) red, (2) 
dark blue, (3) purple, and (4) green are plotted in the ascending order of this parameter. 




Figure 2: Coefficient of restitution as a function of the relative incident velocity: (a) QLV 
Maxwell model, (3 = BvoTR/ho; (b) QLV Kelvin-Voigt model, f = Bvorf/ho. 

Fig. 1^ shows the behavior of the coefficient of restitution as a function of the relative 
incident velocity for (a) QLV Maxwell and (b) Kelvin-Voigt models, respectively. It is to 
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note that with the increase of the relative incident velocity, the coefficient of restitution 
increases for the QLV Maxwell model and decreases for the QLV Kelvin-Voigt model. 
Recall that for the linear Maxwell and Kelvin-Voigt models, the coefficient of restitution 
is independent of the incident velocity. 




Figure 3: Coefficient of restitution as a function of the reciprocal relative mass of the 
impactor: (a) QLV Maxwell model, a = AEqt^/( mho)-, (b) QLV Kelvin-Voigt model, 
a' = AEoot'^/( mho). 


The behavior of the restitution coefficient as a function of the reciprocal relative im¬ 
pactor mass for (a) QLV Maxwell and (b) QLV Kelvin-Voigt models is presented in Fig.|^ 
It is readily seen that with the increase of the relative mass of the impactor, the coefficient 
of restitution increases for the QLV Maxwell model and decreases for the QLV Kelvin- 
Voigt model. Observe that this general tendency for the coefficient of restitution is in 
complete agreement with the behavior of the restitution coefficient in the limiting cases, 
= 0 and {3' = 0, as it follows from formulas (42) and (55). 

In view of the opposite behavior of the the restitution coefficient in the QLV Maxwell 
and Kelvin-Voigt models, it makes sense to consider the behavior of the duration of 
the impact process when the contact between the impact and the tissue specimen takes 
place. Fig. shows the behavior of the relative contact duration as function of the 
relative incident velocity for (a) QLV Maxwell and (b) QLV Kelvin-Voigt models. The 
dependence of the relative contact duration on the reciprocal relative impactor mass is 
plotted in Fig. It is of interest to observe that the relative contact duration for the both 
QLV models decreases with the increase of either relative incident velocity or reciprocal 
relative impactor mass. 

Further, the QLV Zener impact model contains an additional dimensionless parameter 
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Figure 4: Relative contact duration as a function of the relative incident velocity: (a) QLV 
Maxwell model, Tc = tc/rR and {3 = BvoTR/ho] (b) QLV Kelvin-Voigt model, Tc = tc/r'R 
and f3' = Bv^r'i^/ho. 




Figure 5: Relative contact duration as a function of the reciprocal relative mass of the 
impactor: (a) QLV Maxwell model, Tc = t^TR and a = AEQT^/{mho) (b) QLV Kelvin- 
Voigt model, Tc = tc/'TR and a' = AE^t'^/{ mho). 

p. The behavior of the main impact parameters for the QLV standard solid model is 
studied for the following different values of the elastic moduli ratio p\ 0.1, 0.3, 0.6, 0.9, 
and 0.0, corresponding to the limiting case of QLV Maxwell model. In Fig. [6l-[TH in the 
case (a) we use the range of reciprocal relative impactor mass a from 0 to 5 with the 
hxed relative incident velocity (3 = 0.4, while in the case (b) the range of relative incident 
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Figure 6: Coefficient of restitution in the QLV standard solid model as a function of: (a) 
the reciprocal relative mass of the impactor a with the fixed relative incident velocity 
/3 = 0.4 and (b) the relative incident velocity (3 with the fixed reciprocal relative mass of 
the impactor a = 0.4 for different values of the elastic moduli ratio p. 


velocity (3 from 0 to 5 is employed with the fixed reciprocal relative impactor mass a = 0.4. 

As a result of numerical calculations, we observe a nonmonotonic behavior of the 
coefficient of restitution for the QLV Zener model. In Fig. (a) and (b), the initial 
dropping part of each curve (except the limiting case of QLV Maxwell model) is followed 
by the ascending part of the curve after passing through the minimum point, which shifts 
to the right with increasing p. 

The nonmonotonic of the coefficient of restitution is reflected by a non-trivial depen¬ 
dence of the relative contact duration on the parameter p (see Fig. where the curves 
cross over). Nevertheless, the relative contact duration decreases with the increase of 
either relative incident velocity or reciprocal relative impactor mass (see Fig. 7). 

The behavior of the relative time to displacement peak value (see Fig. 8) and the 
relative peak value of impactor displacement (see Fig. is similar and, in particular, 
the two quantities both decrease with increasing the elastic moduli ratio (the direction 
of the increasing p is shown by the wide arrow). Note that the dimensional peak value 
of impactor displacement is given by Xm = VQTR^rni.'Tm)i therefore will increase with 
increasing uq- 

Fig. IT shows that the time to contact force peak value increases with the increase 
of the elastic moduli ratio. At the same time, the peak value of contact force in the 
QLV standard solid model monotonically increases when either the impact mass m or the 
incident velocity Uq increases (see Fig. 11 (a) and Fig. 11 (b), respectively). 
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Figure 7: Relative contact duration Tc = in the QLV standard solid model as a 

function of: (a) the reciprocal relative mass of the impactor a with the fixed relative 
incident velocity (3 = 0.4 and (b) the relative incident velocity (3 with the fixed reciprocal 
relative mass of the impactor a = 0.4 for different values of the elastic moduli ratio p. 




Figure 8: Time to displacement peak value = tm/Tfi in the QLV standard solid model 
as a function of: (a) the reciprocal relative mass of the impactor a with the fixed relative 
incident velocity f3 = 0.4 and (b) the relative incident velocity f3 with the hxed reciprocal 
relative mass of the impactor a = 0.4 for different values of the elastic moduli ratio p. 
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Figure 9: Peak value of impactor displacement imiTm) = Xm/ivoTn) in the QLV standard 
solid model as a function of: (a) the reciprocal relative mass of the impactor a with the 
hxed relative incident velocity /? = 0.4 and (b) the relative incident velocity f3 with the 
hxed reciprocal relative mass of the impactor a = 0.4 for different values of the elastic 
moduli ratio p. 




Figure 10: Time to contact force peak value tm = tulTR in the QLV standard solid model 
as a function of: (a) the reciprocal relative mass of the impactor a with the hxed relative 
incident velocity (3 = 0.4 and (b) the relative incident velocity (3 with the hxed reciprocal 
relative mass of the impactor a = 0.4 for diherent values of the elastic moduli ratio p. 


5 Discussion and conclusion 


First of all, let us compare the predictions of the two limit QLV models for same specihc 
values of the dimensionless parameters. It is of interest that the two models exhibit dif- 
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Figure 11: Peak value of contact force = FiMho/^AEo) in the QLV standard solid 

model as a function of: (a) the reciprocal relative mass of the impactor a with the hxed 
relative incident velocity (3 = 0.4 and (b) the relative incident velocity f3 with the hxed 
reciprocal relative mass of the impactor a = 0.4 for different values of the elastic moduli 
ratio p. 


ferent trends in regard to the coefficient of restitution. Namely, for the QLV Maxwell and 
Kelvin-Voigt models, the quantity e* increases and decreases, respectively, with increasing 
the impactor velocity (see Fig. [^. At the same time, the dependences of the coefficient of 
restitution on the relative impactor mass are similar to those of the corresponding linear 
models (see Fig. [^. 


It is clear that the nonlinear description of the elastic stress function proposed by Fung 


(1981) plays an important role in the qualitative behavior of the solution of the impact 
problem in the QLV Zener model. In particular, introducing the non-linearity parameter 
B, we make the behavior of the coefficient of restitution to be dependent on the impactor 
speed at incidence. Moreover, the choice of the stress-relaxation functions contributes to 
the model’s qualitative performance as well. Observe that, because both parameters f3 
and f3' are proportional to Vq and i?, the effect of increasing non-linearity parameter B 
on the coefficient of restitution e* will be the same as the effect of increasing uq. 

In the case under consideration, the coefficient restitution for the QLV Zener model 
in the range of small values of the incident velocity vq or the reciprocal relative impactor 
mass 1/m slightly decreases similar to the tendency observed in the case of QLV Kelvin- 
Voigt model, while for other considered values of the parameters (3 and a the coefficient 
restitution for the QLV Zener model increases similar to the QLV Maxwell model. It 
can be suggested that the crossing curves of contact duration observed in Fig. [^for the 


17 
































increasing elastic moduli ratio p in the QLV standard solid model can shed a light on the 
behavioral transition from one pattern to the other, corresponding to the limiting cases. 
However, it should be taken into account that while the QLV Maxwell model is recovered 
from the QLV Zener model in the limit as p —)■ 0, the passage from the QLV Zener model 
to the QLV Kelvin-Voigt model requires a certain renormalization as it was shown in 
Section 13.21 

It was shown (Argatov and Mishuris, 2011, Argatov, 2013) that the short-time re¬ 
sponse of a linear biphasic layer, whose deformation is described in the framework of 
the asymptotic model given by Ateshian et ah (1994), in blunt impact is mathematically 
(approximately) equivalent to that of the Maxwell model. Hence, one can anticipate 
that applying the QLV Maxwell model to articular cartilage is more appropriate than the 
choice of of the QLV Kelvin-Voigt model. However, the validity of the QLV models for 
articular cartilage under physiological load magnitudes should be verihed experimentally 
at least for some range of the model parameters. 

Finally, we should point out that there are strong limitations of the assumption of 
quasilinear elasticity for soft biological tissues. In particular, it has been shown that when 
htting the quasilinear viscoelastic law to relaxation curves, creep curves are reasonably 


well reproduced, while hysteresis under cyclic loading is grossly underestimated (Bischoff 


et ah, 2004). Note also (Provenzano et ah, 2002) that one of the known major drawbacks 


of Fung’s QLV model is that its time dependence is independent of strain, therefore 
predicting the same relaxation rate regardless of applied strain. 
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